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We study coherent control in the vicinity of metallic nanostructures. Unlike in the case of control 
in gas or liquid phase, the collective response of electrons in a metallic nanostructure can significantly 
enhance different frequency components of the control field. This enhancement strongly depends 
on the geometry of the nanostructure and can substantially modify the temporal profile of the local 
control field. The changes in the amplitude and phase of the control field near the nanostructure are 
studied using linear response theory. The inverse problem of finding the external electromagnetic 
field to generate the desired local control field is considered and solved. 



Last two decades have witnessed the birth and rapid 
progress in manipulation of quantum systems, such as 
atoms and molecules 0, 0, 0], semiconductor quantum 
dots [i| or even complex biological systems @ by means 
of optimally shaped ultrashort laser pulses. Coherent 
control is usually achieved in homogeneous gas or liquid 
phase by the direct coupling of electronic or molecular de- 
grees of freedom of the system being controlled to the ex- 
ternal electromagnetic fields. However, it is known that 
most of the important chemical reactions in nature and 
technology take place at surfaces or interfaces. For exam- 
ple, at the surface of a catalyst (usually metal) breaking 
and formation of chemical bonds takes place on a fem- 
tosecond time scale and observation of the dynamics can 
be performed, e.g., using a pump-probe technique Q. 
The induction and control of chemical reactions, such as 
oxidization of CO on a metal surface, using ultrashort 
laser pulses is very important from scientific and techno- 
logical points of view 0, 0] • 

For coherent control of atoms or molecules near metal- 
lic surfaces, metallic clusters or other nanostructures 
the situation becomes very different from the control in 
homogeneous dielectric media. The interaction of va- 
lence electrons with external field may lead to collective, 
geometry-dependent excitation modes and produce lo- 
cal fields few orders of magnitude larger than the inci- 
dent control field 0|. In this case the direct coupling of 
the control field to the controlled molecule may be much 
weaker than indirect coupling via interaction with the ex- 
cited plasmon modes. Note, that one has to distinguish 
between the static spatial localization of the induced field 
near sharp edges of metallic nanostructure and excitation 
of spatially localized plasmonic modes. The interaction of 
electromagnetic field with a metallic nanostructure, un- 
like in the case of homogeneous dielectric medium, cre- 
ates highly inhomogeneous induced local field distribu- 
tion. The spatial variations of the field intensity may 
be significant on a scale as small as one nanometer |9j. 
The effect of the presence of a metallic nanostructure is 



not limited to the local field enhancement only, but may 
also result in the broadening of a state linewidth in a 
quantum system due to dynamical screening effects. For 
example, recent measurements demonstrate [Toj j reduc- 
tion of fluorescence lifetime of dye molecules (from 2.8 
ns to less then 1 ns) in the presence of nanostructured 
metal surface. 

Since ab intio simulations of a controlled quantum sys- 
tem and electrons in the nanostructure may be computa- 
tionally extremely expensive, most of the previous mod- 
eling of temporal characteristics of the local fields and 
associated quantum dynamics was typically performed 
for oversimplified models, with the focus on some partic- 
ular aspects of the problem. For example, a quantum me- 
chanical treatment of time-dependent screening and asso- 
ciated transient effects due to sudden creation of a charge 
in a close vici nity of two dimensional electron gas, was 
considered in llj . Other simulations, involving more re- 
alistic geometries than an idealized two-dimensional elec- 
tron gas, are usually performed using classical Finite Dif- 
ference Time Domain (FDTD) approach. One can men- 
tion coherent control of nanoscale localization of optical 
excitations in nanosystems [9(, where classical equations 
with local phenomenological dielectric constant e(r) were 
used. However, the validity of local classical response 
theory is questionable for relatively small or highly inho- 
mogeneous nanostructures, where quantum mechanical 
effects may play a critical role. In the quantum limit one 
needs to take into account the discreteness of the energy 
spectrum, non-local electron density response, dynamical 
screening effects and tunneling [l2j . 

As can be seen, the coherent control of atoms and 
molecules localized near metallic nanostructures is much 
more complex than in a gas or liquid phase. It is impor- 
tant for the theory to capture the fact that the incident 
electromagnetic field is coupled both to single particle 
and collective plasma excitations in the nanostructure, 
and the field in the nanostructure is neither completely 
screened, nor purely transverse. To develop the theory 
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of the coherent control at the nano (i.e. subwavelength) 
scale, we need to map out spatio-temporal or spatio- 
spectral amplitude as well as the phase of the local field 
in relation to the external spatially-uniform, but tempo- 
rally shaped excitation field. More importantly, the in- 
verse problem of finding the external field corresponding 
to the required local control field for the best performance 
needs to be addressed. It may seem that the modification 
of the control field in the vicinity of a metal nanostruc- 
ture makes coherent control very difficult or impossible. 
However, we will show that one can successfully generate 
local control fields with the prescribed temporal behavior 
and take advantage of local field enhancement. 

In this study we consider the problem of coherent con- 
trol in the presence of a doped semiconductor nanostruc- 
ture and propose to use the local field enhancement due 
to geometrical and resonant plasmon excitation for op- 
timal control of atoms and molecules. Relatively long 
relaxation and dephasing times were reported for heavily 
doped semiconductor heterostructures [l3j ]. The carrier 
concentration in the nanostructure is chosen to generate 
plasmon resonance around 0.8 eV, which corresponds to 
the wavelength A = 1550 nm. The experimental tech- 
niques for generation, shaping and detection of ultrashort 
pulses at this wavelength are readily available. Also, this 
wavelength range is important for excitation and optimal 
control of molecular wave packet formation using a two- 
color pump-probe scheme 
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The Fermi wavelength 
at this carrier concentration is w 10 nm, comparable to 
the size of the nanostructure (several tens of nanometers- 
within capabilities of modern nanofabrication methods), 
which motivates us to use quantum nonlocal linear re- 
sponse theory Using the theory [jjj], we examine 
the inhomogeneous distribution of the induced field, and 
study how the control pulse intensity and phase are trans- 
formed due to the plasmonic response of a small custom- 
engineered nanostructure. 

In the self-consistent treatment the induced charge 
density in the nanostructure Pi n d( r )£) is given by: 



dt' / dT , x (T,T , ,t,t')4 >tot {r',t'), (1) 



Pind(r,£) 



where </>tot(r, t) = </> xt(r, t) + </>i n d(r, £) is the total poten- 
tial, and X o(r,r',M') = -6{t-t') < p(r,t), p(r' ,£') > 
is the electron polarizability, and the statistical average 
is performed over the equilibrium state. The function 
8(t — t') ensures the causality of the electron response. 

We assume the density-density response function is 
translation invariant in time variables Xo{ r i r 'i t> = 
Xo(r, r', £ — £'), but not in the space variables. By making 
Fourier transform in time domain and assuming Random 
Phase Approximation (RPA), one can express Xo( r , r ' , w ) 
in terms of the eigenenergies and eigenfunctions of the the 



unperturbed Hamiltonian: 

X0 (r,r» = X; J^li^} ^ *(r)^(rOV>*Wi(r), 



(2) 

where / is the Fermi distribution function, and 7 ac- 
counts for the finite width of the quantum levels in the 
nanostructure. The eigenproblem for electrons confined 
in the nanostructure Hipi(r) = Eitpi(r) is solved using 
numerical diagonalization. 

The induced potential 0i n d(r,u;) is determined from 
the self-consistent integral equation [l§] 



</>ind(l", Uj) 



Xo (r', r" , W )&at (r" , w) V c ( | r - r' | )dr' dr" ,(3) 



where Vc(|r — r'|) is the Coulomb potential. Here we 
make use of the fact that the nanostructure is much 
smaller than the wavelength of the incoming field, so 
the retardation effects can be neglected. In order to find 
</>ind(r, uS) we discretize Eq.Q on a real-space cubic mesh 
and solve it numerically. The equation Eq. ([3]) establishes 
one-to-one correspondence between the homogeneous ex- 
ternal control field E oxt (o;) and the total (local) field 
Etot(r,o>) at a given point r. Note that this correspon- 
dence is established using statistical average and has no 
uncertainties related to quantum dynamics of electrons 
for nanostructures with large number of electrons. 

Since we work in a weak-field limit and assume linear 
response of the system, there is no mixing or creation 
of new frequency components, which were not originally 
present in the incoming control field. This linear response 
approximation significantly simplifies the analysis in that 
the response of the system can be computed for each 
frequency separately in the spectral domain: 



E to t(ro,w) = Z(r , w)E 0Xt (u), 



(4) 



where i"o is the position of the controlled molecule, and 
for a small metallic nanostructures the complex valued lo- 
cal response kernel Z(ro,w) is non-zero for any to. Thus, 
the expression Eq.Q can be inverted, and the optimal 
external field can be computed given the required local 
control field E to t(ro, oj). 

To illustrate our approach we consider a finite nanos- 
tructure with dimensions L x L x <i, L = 33 nm, d = 
10 nm, made of a doped semiconductor with the car- 
rier concentration n = 9.5 x 10 19 cm~ 3 with a bow-tie 
shaped hole, FigQTa). Nanostructures of such shape are 
routinely manufactured for local field enhancement [HI]. 
There are total 1042 electrons confined in the nanostruc- 
ture. The electron effective mass is taken to be 0.2m e , 
where m e is the bare electron mass [la ]. The carrier 
concentration is chosen to generate plasmon resonance 
around 0.8 eV. 

In our simulations we set the phenomenological damp- 
ing parameter to 7 = 3.8 meV 13j. We performed 
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FIG. 1: (a) Normalized local field intensity Etot(i")| 2 /|E 

ext | in 

the nanostructure with a bow-tie shaped hole, for static (u> = 0) 
external field. Note strong field enhancement near the sharp edges, 
(b) Normalized local field intensity |Etot(i"o)| 2 /|E ox t| 2 (solid line) 
and phase (dashed line) in the center of the nanostructure as a 
function of photon energy of the external field hv. Note many local 
extrema due to multiple geometric resonances in the nanostructure. 
The resonance peak occurs at hv 0.8 eV as the phase reaches 7r/2 
value. 



simulations assuming room temperature T = 300 K ks 
25.8 meV and linearly polarized external field with the 
polarization direction parallel to the x axis. 

Figure [TJa) shows spatial distribution of the normal- 
ized total field intensity |Etot(r)| 2 /|E cx t| 2 in the static 
(u> — > 0) limit. Note strong field enhancement in the 
vicinity of the bow-tie sharp edges. We assume that con- 
trolled molecule is localized in the center of the hole, in 
the point of local field maximum, which is our target 
region. 

We performed calculations of the local field E to t in the 
target region for different frequencies of the external field, 
Fig. QJb) . Note that the main plasmonic resonance peak 
occurs at u ~ 0.8 eV as the phase reaches ir/2 value. 
The frequency dependence has many local extrema due 
to multiple geometric resonances in the nanostructure. 
Using the spectral response data of Fig. [Hb) i the tem- 
poral structure of the local field induced by an arbitrary 
incident field can be uncovered, Eq. (j4|) , by Fourier trans- 
forming the product of the response of the nanostruc- 
ture and the incident field complex spectrum. Similarly, 
the external control field in the frequency domain can be 
computed simply by dividing the given local field com- 
plex spectrum by the nanostructure response function. 

Some examples are shown in Fig. [21 A 10 fs transform- 
limited Gaussian pulse, Fig.[5Ja,b) incident on the nanos- 
tructure will result in a local field shown in Fig. [2Jc,d). 



Note the substantial difference in the shape of the lo- 
cal field as compared to the incident field, particularly 
the parasitic "ringing" following the main pulse. This is 
not surprising, since in the simplest view the collective 
response of the electrons in the nanostructure comprise 
a damped oscillator with one or more eigenmodes, the 
frequencies and amplitudes of which are determined by 
the geometry and electron concentration. The solution 
to the inverse problem of obtaining the local field of the 
form shown in Fig. [2[a,b) is displayed in Fig. [2je,f) and 
is easy to understand by noticing that the spectral phase 
of the incident field is simply the inverted phase of the 
nanostructure response. Similarly the spectral amplitude 
structure of the incident field compensates for the spec- 
tral variations of the nanostructure response to result in 
a smooth amplitude spectrum of the local field. A con- 
sequence of this is that the "ringing" on the drive pulse 
is exactly out of phase with the nanostructure "ringing" , 
so that the two contributions cancel each other. 




FIG. 2: (a) Time dependence of a Gaussian external field pulse 
with duration 10 fs. (b) Spectral intensity and phase of the in- 
coming pulse shown in (a). Note, the phase of the pulse is set to 
zero, (c) Time dependence of the local field excited in the target 
region by the pulse (a), (d) Spectral intensity and phase of the 
local field shown in (c). (f) Time dependence of the external field 
that excites the total field in the target region, same as in (a), (e) 
Spectral intensity and phase of the incoming pulse shown in (f). 

In conclusion, we have demonstrated theoretically 
the feasibility of optimal coherent control near metallic 
nanostructures, and how to use the resonant properties 
of the nanostructure to generate strongly enhanced, spa- 
tially and temporally localized control field. We have 
shown that the local field shape deformation and para- 
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sitic "ringing" can be successfully eliminated. For rela- 
tively weak fields one can readily find an external control 
field, that for the particular geometry and given carrier 
concentration in the nanostructure will generate the nec- 
essary local control field at a given point. The solution 
of the inverse problem always exists and does not depend 
on the shape of the optimal local pulse. In practice, the 
optimal control pulse can be found using standard closed 
loop techniques [l7|, because the exact response of each 



given nanostructure will be nearly impossible to mea- 
sure. Moreover, we studied the frequency dependence of 
the local field at various points within the nanostructure. 
Although the inhomogeneous local field amplitude may 
vary significantly from point to point, we discovered that 
the phase of the field is much less sensitive to spatial vari- 
ations. This provides additional tolerance to the position 
fluctuations of the controlled molecule. 

In our simulations we used a model of non-interacting 
electrons trapped in the nanostructure. Such simplifi- 
cations will potentially lead to underestimation of the 
Landau damping [III, resulting in narrower resonance 
peaks. However, the position of the main resonance 
peaks obtained, for example, using the ab initio Den- 
sity Functional Theory, should be similar to our simpli- 
fied model [18[. The applicability of the linear response 
theory, used in this study, is limited to relatively weak 
fields. For many optimal control applications this as- 
sumption holds. However, in some cases, the field inten- 
sity can be high enough, resulting in significant tempera- 
ture rise and strong nonequilibrium electron distribution 
which will require the use of the nonlinear response 
theory. Further, the breakdown of the linear response 
approximation will make the response kernel Z(ro, to) de- 
pendent on the amplitude of the external field. In this 
case one cannot guarantee the inversion of Eq. (j4]) for an 
arbitrary optimal local field. 

The electron-electron, electron-phonon scattering and 
other processes, having unpredictable, stochastic nature 
lead to thermalization of the electron gas in the nanos- 
tructure and eventually to the loss of coherence in the 
system [l3j]. However, if the controlled molecule is not in 
the direct contact with the nanostructure, these processes 
would not limit the possibility of the optimal control. Be- 
cause of the relatively large number of electrons in the 
nanostructure, the local electric field has no uncertain- 
ties related to the quantum dynamics of electrons, and 
the nanostructure effectively acts as a linear filter. In 
this case the coherent control is limited by the intrinsic 
decoherence times of the controlled molecule. However, 
as we mentioned above, these times may be decreased 
due to the dynamical screening effects [1 01 ] . 

Since the temporal profile of the local control field 
keenly depends on the position of the controlled molecule 
on the nanostructure, one needs to address possible 
mechanisms and controllability of spatial localization of 
the controlled molecules. It was theoretically demon- 



strated that the plasmon-generated spatially inhomoge- 
neous field between metal nanoparticles can serve as a 
trap for molecules [n| without direct contact to the 
nanoparticles, which otherwise would lead to strong deco- 
herence. It was also shown 12, 19[ that the field enhance- 
ment is maximum in the space between the nanoparticles. 
Thus, properly arranged arrays of metal nanoparticles 
may offer both trapping and field enhancement capabili- 
ties. 

This work was performed, in part, at the Center for 
Integrated Nanotechnologies, a U.S. Department of En- 
ergy, Office of Basic Energy Sciences user facility. Los 
Alamos National Laboratory, an affirmative action equal 
opportunity employer, is operated by Los Alamos Na- 
tional Security, LLC, for the National Nuclear Security 
Administration of the U.S. Department of Energy under 
contract DE-AC52-06NA25396. 
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